Formation of Heavy Meson Bound States by Two Nucleon Pick-up Reactions 
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We develop a model to evaluate the formation rate of the heavy mesic nuclei in the two nucleon 
pick-up reactions, and apply it to the ®Li target cases for the formation of heavy meson— a bound 
states, as examples. The existence of the quasi-deuteron in the target nucleus is assumed in this 
model. It is found that the mesic nuclei formation in the recoilless kinematics is possible even for 
heavier mesons than nucleon in the two nucleon pick-up reactions. We find the formation rate of 
the meson— a bound states can be around half of the elementary cross sections at the recoilless 
kinematics with small distortions. 
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I. INTRODUCTION 

Meson-Nucleus systems are one of the most interest- 
ing laboratories to study the meson properties at finite 
density and to explore the symmetry breaking pattern 
of QCD and its partial restoration in nucleus Es- 
pecially, the studies of the bound states of meson and 
nucleus have the following advantages; (i) the selective 
observation of the meson properties is possible by making 
use of the fixed quantum numbers of the bound states, (ii) 
the meson properties inside the nucleus can be observed 
clearly only with relatively small contamination from the 
vacuum processes, and (iii) the system is quasi-static and 
the time dependent dynamical evolution of the system is 
irrelevant. These features are different from other meth- 
ods based on the scattering and collision processes. On 
the other hand, the information obtained from the ob- 
servation of the bound states are limited in p < po and 
T = region in the QCD phase diagram. 

We have studied so far the physical interests, the struc- 
tures, and the formation reactions of the various kinds of 
the meson- nucleus systems 0-11]. Within these studies, 
the most exciting and successful results were obtained 
by the observation of deeply bound pionic atoms in the 
one nucleon pick-up (c?,'^He) reactions 0, 0, 0| ■ We have 
also considered other one nucleon pick-up reactions like 
{l:P) 9 and (tt, A^) Q for the mesic nuclei formations. 
The one nucleon pick-up reactions are found to be useful 
for the mesic nucleus formation, however, they require 
the large momentum transfer for heavy meson produc- 
tion which is one of the main obstacle to observe bound 
states. Thus, we think that we need to develop new 
methods for bound state formations to extend our stud- 
ies to other meson-nucleus systems, especially for heavier 
mesons. Actually, we are interested in the heavy meson- 
nucleus systems such as 7/(958) meson for the studies of 
C/a(1) anomaly effects [ICHL^ I. <j) meson for ss compo- 
nents of nucleon and OZI rule at finite density [isj . and 
D mesons for the charm meson properties in nucleus [3] . 
Thus, we would like to study the two nucleon pick-up re- 



actions theoretically as a possible method suited to form 
the heavy meson nucleus bound systems. 

One of the most serious problems in the formation re- 
actions of the heavy meson nucleus systems are the large 
momentum transfer as mentioned above. It is known 
that the matching condition of momentum and angular 
momentum transfers plays an important role to deter- 
mine the largely populated subcomponents and it is also 
known that the best choice for our purpose is the to- 
tal angular momentum transfer J = state formation 
in the recoilless kinematics in many cases. In the one 
nucleon pick-up reactions which we have mainly consid- 
ered so far, the recoilless kinematics can not be satisfied 
for the formation of heavier meson bound states than 
the nucleon mass because of the large mass of the me- 
son. Thus, we consider the two nucleon pick-up reac- 
tions in this article to investigate the possibility to extend 
our study to heavier meson region using these reactions. 
Actually, there was an attempt to observe the ?7-mesic 
state in the two nucleon pick-up ^^Al(p,'^He) reaction at 
COSY-GEM M- 



II. EFFECTIVE NUMBER FORMALISM FOR 
THE QUASI-DEUTERON IN NUCLEUS 



We formulate the formation cross section of the heavy 
meson bound states by the two nucleon pick-up reactions. 
As we will see below, our model is so simple that it can be 
applied generally to the heavy meson bound state forma- 
tions by the two nucleon pick-up reactions such as (7, d) 
and (p,^He). 

We apply the effective number approach, which has 
been used for the studies of the meson-nucleus bound 
states Q, to evaluate the formation rate of the bound 
systems in the two nucleon pick-up reactions. In the 
effective number approach, the formation cross section 
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by the two nucleon pick-up reactions can be written as, 
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where {dujd^Y^^ is the elementary cross section of the 
meson production, and F the width of the meson bound 
states. The all combinations of the final states, labeled by 
/, are summed up to evaluate the inclusive cross section. 
The energy transfer Ai? of the reaction in the laboratory 
frame is defined as, 

l\E ^ {Tf + Mf) + {M -B)- {T, + M,) - {2Mn - S2n), 

(2) 

where Tf and Mf are the kinetic energy and the mass of 
the emitted particle, Ti and Mi the kinetic energy and 
the mass of the incident particle, and M the mass of 
the produced meson. The meson binding energy B and 
the two nucleon separation energy S'2Ar from the target 
nucleus are determined for each bound level of meson and 
excited level of the daughter nucleus. Here, we neglect 
the recoil energy of the nucleus. 

The momentum transfer q of the reaction is defined as. 



q = Pi -Pf, 



(3) 



and shown in Fig. [T] for the (f> meson formation case, as 
an example, in the heavy target for the (7, d) reaction to- 
gether with that for one-nucleon pick-up (7, N) reaction 
as functions of the momentum of the incident photon. 
We also show the momentum transfers for the (p,'^He) 
and {p, d) reactions for the (j) meson formation in Fig. [2] 
for comparison. As we can see from the figures, the two- 
nucleon pick-up reactions satisfy the recoilless condition 
q = at a finite incident momentum, while large momen- 
tum transfer is unavoidable for the one-nucleon pick-up 
reactions. We show in Fig. [3] the incident beam momenta 
p for the four reactions required to produce the meson 
with an effective mass M* in recoilless kinematics. The 
effective mass M* is defined as M* = M — B with the 
meson mass M and the binding energy B. We find clearly 
that the meson with larger effective mass than nucleon 
can not be produced in recoilless kinematics by the one 
nucleon pick-up reactions. In two nucleon pick-up reac- 
tions like (7,c?) and (p,'^He), on the other hand, we can 
produce the heavier meson states like 77'(958), (/)(1020), 
and ai(1260) meson bound states in the recoilless kine- 
matics. 

The theoretical calculation of the two nucleon pick-up 
reactions is rather difficult in general. In addition, we ex- 
pect complicated nuclear excitations in daughter nuclei 
with two nucleon holes, which will prevent us from the 
clear identification of meson bound states. Here, we con- 
sider the specific nucleus ^Li as the target which has the 
well-developed cluster structure of a -I- d in the ground 
state. The probability of the a + d component of the 
ground state of ^Li is reported to be 0.616 in Ref. (T6j 
and 0.73 in Ref. 17|. In the reaction considered in this 
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FIG. 1. Momentum transfer q of the forward (7,d) reaction 
for the formation of (ji(1020) meson bound states in the heavy 
target nucleus plotted as functions of the incident photon mo- 
mentum for four values of the gap (B^ — Sd) between the <j!> 
meson binding energy and the deuteron separation energy 
Sd as indicated in the figure. Momentum transfer of the one- 
nucleon pick-up (7, A'^) reaction is also shown for comparison. 



800 



600 



400 



200 




3 4 
Pp [GeV/c] 

FIG. 2. Same as Fig. [T]except for the proton induced (p,'^He) 
and (p, d) reactions. 



article, the momentum transfer to the initial deuteron 
wavefunction in ®Li and the final meson wavefunction in 
the mesic nucleus is considered to be small near the re- 
coilless kinematics and the quasi-deuteron picture in the 
^Li target is expected to be a good approximation. Thus, 
by considering the nuclei like ^Li which have the large 
component of quasi-deuteron as targets, we can evaluate 
the reaction rate in a simple way and expect to have the 
simple structure of the formation spectra of mesic nu- 
clei in the two nucleon pick-up reactions. In our model 
considered, we treat ^Li as the bound state of the alpha 
particle and the deuteron to evaluate the reaction rate. 
Thus, B in Eq. ^ is the binding energy of a heavy me- 
son and a particle in the final state and S2N is fixed to 
be S2N = 1.47 [MeV] from the mass gap of initial and 
final nuclei as S2N = {Ma + Mj) — Men. The elemen- 
tary cross section [dcr/dVlY^'' for the meson production 
appearing in Eq. ([1]) is that of the i -\- d ^ f + meson 



3 




'0 400 800 I I 1200 1 1600 

Mri Mr^'M^ Mai M* [McV] 




r[fm] 



FIG. 3. Incident beam momenta p required to produce meson 
in recoilless kinematics in one nucleon pick-up ((7, A^), {p, d)) 
and two nucleon pick-up ((7, d), (p,"^IIe)) reactions are plotted 
as functions of meson effective mass M* which is defined as 
M* = M — B with the meson mass M and the binding energy 
B. The in- vacuum masses of rj, rj' , (p and ai are indicated in 
the figure by arrows. 



reaction with the incident particle i and the emitted par- 
ticle /, which should be evaluated from experiments as 
previous cases Here, the emitted particle / will be 

^He for the proton induced (i = proton) case, and d for 
the 7 induced {i — photon) case, respectively. 

In our model considering the ^Li nucleus as the bound 
state of alpha particle and deuteron, the effective number 
A^eff of the ^Li(i, /)aCS)meson reaction can be written as. 
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where (f>i^ (r) is the wavefunction of the meson bound 
state and ipiaif) that of the deuteron bound to a in the 
^Li target. Xii^) ^-nd X/C*") ^-re the incident and the 
emitted particle wave functions in the scattering states, 
respectively. We assume plane waves for Xi and X/ in 
this article. The distortion effects to Xi ^-nd Xf depend 
on the incident and emitted particles however they 
are known to be relatively small for the cases satisfying 
the matching condition [J]. The deuteron wavefunction 
ipi^ in ^Li is determined to reproduce the momentum 
distribution reported in Ref. based on the analysis 
of the ^Li(e, e'd)^He reaction. We calculate the tpi^ by 
solving the Schrodinger equation with the Woods-Saxon 
type potential. 



Uir) 



Uo 



l + cxp[(r-i?)/a]' 



(5) 



and adjust the potential depth Uq and the radius pa- 
rameter R to reproduce the momentum distribution p(jj) 
reported in Ref. [13] ■ p{p) is defined as. 
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FIG. 4. Calculated density distribution of the radial part 
Ri^(r) of the relative wavefunction tjji^{r) of a and deuteron 
in ®Li nucleus. 



where |V'/d('")P normalized to be 1 in the coordinate 
space as usual. The potential parameters used here are 
fixed to be i? = 2.0 [fm], a = 0.5 [fm], and Uq = -75 
[MeV] . The calculated wavefunction is shown in Fig |3] 
and momentum distribution in Fig. [5] The wavefunction 
in Fig. H] co rresponds to the 2s bound state as indicated 
in Ref. [13] because of the Pauli effect to nucleons which 
forbids the alpha and deuteron clusters to be in the rel- 
ative Is state. The calculated momentum distribution 
reproduce the PWIA result in Ref. flf\ reasonably well 
as shown in the Fig. [SJ In Fig. [5l we multiply a factor 
0.73 to our results to correct the overall normalization 
of p{p) to be the same as Re f. IfTII . It should be noted 
that the PWIA results in Ref. [14] ignoring the distortion 
effects are the quantity which should be compared with 
our results calculated by Eq. 

The bound meson wavefunctions (f>i^ (r) in the final 
state are calculated by solving the Klein-Gordon equa- 
tion with the optical potential in the Woods-Saxon form 
as in Eq. ([5]). Wc fix the meson mass to be M = 1020 
[MeV] as (f) meson, which is heavier than nucleon and can- 
not be formed in the recoilless kinematics in one nucleon 
pick-up reactions. Since mesons can be absorbed by the 
nucleus generally, the potential strength Uq in Eq. ^ is 
considered to be a complex number as Uq = (Vq -I- iWo) 
for mesons. In the present calculation, we consider a 
few potential strengths as examples and study how the 
two nucleon pick-up reaction spectra change according 
to the meson-nucleus interaction. We consider first the 
potential strengths based on the experimental data ob- 
tained by E325 in KEK [l^, where the meson mass 
shift is reported as Am{pQ)/m — —3.4 % and the (p 
meson width in nucleus T^{po) — 15 [MeV]. We adopt 
these numbers as the (j) mesic optical potential and fix 
as {Vo,Wo) = (-34.7,-7.5) [MeV]. We mention here 
that the large in-medium width r0(po) — 80 [MeV] of 
(j) meson was also indicated based on another attenua- 
tion experiment by LEPS at SPring-8 [l9j. which corre- 
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TABLE I. Calculated binding energies and widths of the 
meson— alpha bound states in unit of MeV with the poten- 
tial strength {Vo,Wo) = (i) (-34.7,-7.5), (ii) (-250, -5), 
and (iii) ( — 250, —20) [MeV]. The meson mass is fixed to be 
1020 MeV. 
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FIG. 5. Momentum distribution of deuteron in "^Li obtained 
by our model (dashed line) and by the analysis of the ex- 
perimental data [13] (solid line). The correction factor 0.73 
is multiplied to the results of our model to satisfy the same 
normalization as in Ref. U3i- 



sponds to the imaginary potential strength Wq = —40 
[MeV]. We have checked numerically that the potential 
with this imaginary strength together with the real part 
strength corresponding to the 3.4 % mass reduction does 
not provide ((> meson bound states in alpha potential. We, 
then, assume a stronger attractive potential for meson- 
nucleus system to estimate the formation rate of heavy 
meson bound systems with strong attractive potential as 
reported in Ref. [ll| for 77'(958) meson in a theoretical 
model. The assumed potential parameters with two dif- 
ferent absorption strength are (Vb, Wq) = (—250, —5) and 
(—250, —20) in unit of MeV. The distribution parameters 
are fixed to be i? = I.ISA^/^ - 0.48 [fm] and a = 0.5 [fm] 
with A = 4 for meson-alpha system. We show in Table U 
the calculated binding energies and widths of the me- 
son bound states for the three different potentials. The 
radial density distributions are also shown in Fig. [6] for 
{Vo,Wo) = (-34.7,-7.5) and (-250,-5) [MeV] cases. 
We found the radial density distributions of Is states of 
the both potentials are much different because of the dif- 
ferent potential strength. 

Since the angular momentum Id of the relative wave- 
function of deuteron and a particle in ^Li is considered 
to be 0, the expression of the effective numbers in Eq. 



state 


(i) 

B.E. 


r 


(ii) 

B.E. 


r 


(iii) 

B.E. 


r 


Is 


0.76 


2.9 


131.5 


8.1 


131.2 


32.5 


2s 






14.5 


2.8 


14.2 


11.2 


2p 






58.0 


5.5 


57.7 


21.8 


3d 






2.9 


2.8 


2.4 


11.3 



1.0 
0.8 



E 

t. 0.6 



0.4 

0.2 



1 1 

_ A 

Is / \ 2p 


1 


1 1 










^^'s 2s 






- ""-i------ 





r [fm] 

FIG. 6. Calculated density distribution of the radial 
part of the wavefunction <;/>;,„ (t") of the meson 

bound state in the a particle for the potential strength 
(Vo,VFo)=(-34.7, -7.5) [MeV] (thick solid line) and (-250, 
-5) [MeV] (thin lines). 



can be simplified as, 

iVeff = E I / e'^'C (OV'o(r-)dr f, (7) 

in the plane wave approximation. We use this expression 
to calculation the effective numbers in this article. 

In the recoilless kinematics, the only s states of the 
meson bound states can be populated because of the or- 
thogonality of the angular part of the wavefunction as 
can be seen in Eq. ([7]). And because of the approximate 
orthogonality of the radial parts of (/)/^ and '00 j the sub- 
stitutional 2s state of the meson is expected be largely 
populated. 



III. NUMERICAL RESULT FOR BOUND STATE 
FORMATION RATE 

To investigate the momentum transfer dependence of 
A'off for each bound state formation, we show in Fig. [7] 
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FIG. 7. Calculated effective numbers are plotted as functions 
of the momentum transfer of the two nucleon pick-up reac- 
tions for ^Li target for the Is meson bound state formation 
with the potential strength {Vo,Wo) = (-34.7,-7.5) [MeV] 
(thick solid line) and the Is, 2s, 2p, and 3d bound states with 
(Vo,Wo) = (-250,-5) [MeV] (thin lines). 



the calculated effective numbers N^s by Eq. ([7]) as func- 
tions of the momentum transfer |q| for two nucleon pick- 
up reactions for ^Li target. Each effective number has 
the characteristic behavior due to the matching condi- 
tion of the momentum transfer and the angular momen- 
tum transfer. As we have mentioned in the previous sec- 
tion, the effective numbers NeS is exactly for 2p and 3d 
states at |q| = because of the orthogonality condition 
of the angular part wavefunction to the s-wave function 
tpQ. And the substitutional 2s state of the bound meson 
have the largest contribution at \q\ = as we expected. 
The contribution of the Is state with weaker real poten- 
tial Vq = —34.7 [MeV] has the stronger dependence on 
q as naturally expected its larger special dimensions as 
shown in Fig. [H] As the momentum transfer increases, 
the 3d bound state formation has the largest contribu- 
tion at 160 < |q| < 270 [MeV/c] and, then the Is bound 
state at 270 < \q\ < 450 [MeV/c]. Overall strength of 
the heavy meson bound state formation cross section be- 
comes smaller for the kinematics with larger momentum 
transfer. 

We show the same effective numbers as functions of the 
incident particle energies for {j,d) and (p, '^He) reaction 
cases in Figs. |8] and HI Because of the different binding 
energies of the meson appearing in Eq. ([5]), the incident 
particle energy which corresponds to recoilless kinematics 
{\q\ = 0) is different for each bound state. We indicate 
the incident particle energy of the recoilless kinematics 
for each state by the arrow in Figs. |5] and [HI The both 
figures show the very similar behavior of Ncg. 

We find that the substitutional 2s state of the meson 
is produced in the almost recoilless condition and has 
the largest contribution to the cross section at pj = 1.6 
[GeV/c] (Tp = 2.7 [GeV]). (In the explanation below. 
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FIG. 8. Same as Fig. [7] except for the plots as functions of 
the incident photon momentum p-y for (7, d) reaction. The ar- 
row indicates the incident photon momentum of the recoilless 
kinematics for each meson bound state formation. 
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FIG. 9. Same as Fig. [7| except for the plots as functions of 
the incident proton kinetic energy for (p, ^He) reaction. 
The arrow indicates the incident proton kinetic energy of the 
recoilless kinematics for each meson bound state formation. 



we indicate the photon momentum p-y in (7, d) reac- 
tion with the corresponding proton kinetic energy Tp in 
the (p,'^He) reaction in the parenthesis.) We mention 
here that the Is state wavefunction with the potential 
strength (V^, VFq) = (-34.7, -7.5) [MeV] has larger spa- 
tial dimension which violates the approximate orthogo- 
nality of the radial part with ■00 a-t q = in Eq. ([7]), and 
thus, the contribution of this Is state also has the large 
contribution at p^ = 1.6 [GeV/c] (Tp = 2.7 [GeV]). At 
p^ = 1.35 [GeV/c] (Tp = 2.1 [GeV]), a little above the 
production threshold of the elementary process, the size 
of the effective number of the 2s and 3d states are sim- 
ilar and larger than those of other states formation. At 
p^ = 1.8 [GeV/c] (Tp = 3.1 [GeV]), the effective num- 
ber of the 2s state formation is still dominant, however 
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FIG. 10. Expected spectra of the forward two nucleon pick-up (7, d) reaction ^Li for the formation of the meson— a bound 
states are plotted as functions of the emitted deuteron kinetic energy for the incident photon momenta p-, — (a) 1.35, (6) 
1.6, (c) 1.8, and (d) 2.0 [GeV/c]. The vertical dashed hne indicates the meson production threshold. Each line is calculated 
with the different optical potential parameters for meson— a system as indicated in (a). The spectrum with the potential 
(Vb, Wo) = (—34.7, 7.5) [MeV] is only shown in (b) and (c) because it is small and invisible in (a) and (d). 



the other contributions of the Is and 3d states forma- 
tion become relatively more important than at p-y = 1.6 
[GeV/c] (Tp = 2.7 [GeV]). The contributions of the Is, 
2s, and 3d state formations are important at Pj = 2.0 
[GeV/c] {Tp = 3.5 [GeV]) and the formation spectrum is 
expected to be a little more complicated than those at 
other photon momenta. We find that N^g for 2s bound 
state formation takes the largest value A^off = 0.606 at 
p^ = 1.59 [GeV/c] (Tp = 2.66 [GeV]) which means that 
the meson bound state formation cross section can be 
about half of the elementary cross section. 

We then calculate the relative strength of the forma- 
tion spectra of the meson bound states in alpha par- 
ticle by the two nucleon pick-up reaction in ^Li target 
at incident photon momenta p-y — 1.35, 1.6, 1.8 and 2.0 
[GcV/c] for the (7,(i) reaction which corresponds to the 
Tp = 2.1,2.7,3.1 and 3.5 [GeV] for the (p,^He) reaction. 
We can expect to observe the different behavior of the 
formation spectra at these energies as expected from the 
energy dependence of the effective numbers. The cal- 



culated results - — -- / (-r=- I for (7,6?) reaction are 
dEdill \d\l J 

shown in Fig. [10] for three different optical potential pa- 
rameters for the meson bound in a particle. 

The expected spectra for the potential strength 
(Vo,Vyo) = (—34.7,-7.5) [MeV] case are simple since 
there is only lightly bound Is state. As shown in Fig. [TU] 
(b) and (c), this Is state is seen as a peak close to this 
meson production threshold for p^ = 1.6 and 1.8 [GeV/c] 
{Tp = 2.7 and 3.1 [GeV]). The contribution of this state 
has so strong q and incident energy dependence as shown 
in Figs. |7][9] and that it becomes smaller than those with 
deeper potential cases and invisible in Fig. [10] (a) and 
(d). 

For deeper potential cases with (Vb, Wq) = (—250, —5) 
and (-250,-20) [MeV], we find that the spectra at 
p^ = 1.6 and 1.8 [GeV/c] (Tp = 2.7 and 3.1 [GeV]) are 
dominated by the 2s state formation and the other con- 
tributions are significantly small. On the other hand, we 
can observe clear peak structures of the Is and Sd states 
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formation in addition to the 2s state at = 1.35 and 
2.0 [GeV/c] {Tp = 2.1 and 3.5 [GeV]). The size of the 

spectra - — — / -— are relatively large at py — 1.6 
dEdill \dil J 

and 1.8 [GeV/c] (Tp = 2.7 and 3.1 [GeV]) where the nio- 
mentum transfer of the two nucleon pick-up reactions is 
small. At = 1.8 and 2.0 [GeV/c] {Tp = 3.1 and 3.5 
[GeV]), the momentum transfer is larger and the size of 
the spectra become smaller rapidly for the larger incident 
momentum and energy. We also show the effects of the 
imaginary part of the optical potential to the formation 
spectra in Fig. (TU] Since we have only one deuteron state 
in the initial nucleus ^Li, the reaction spectra have simple 
structure, especially for {Vo,Wo) = (-34.7,-7.5) [MeV] 
potential case. Thus, the main effects of the absorptive 
potential are found to reduce the peak height with large 
width. The overlap of the resonance peaks due to the 
widths only happens between 2s and 3d states for the 
absorption potential strengths studied here in the two 
nucleon pick-up spectra. 

In case if the imaginary potential is large as Wq — ~40 
[MeV] as indicated based on the data reported for the 
meson in Ref. [19], the height of all peaks in the spectra 
shown in Fig. [TU] becomes lower in inverse proportion to 
the strength of the imaginary potential and the contribu- 
tions of 2s and 3d states for Vq — —250 [MeV] case can 
not be distinguished because of the large widths. 

IV. SUMMARY 

We have considered the two nucleon pick-up reactions 
in this article to investigate the feasibility of the reactions 
of this type to extend our research field of the meson 
nucleus bound systems to the heavier meson region. We 
have developed a model and used the effective number 
approach to evaluate the formation rate. As an example, 
we have applied the model to the heavy meson bound 
state formation in alpha particle with ^Li target. As 
shown in the numerical results, we have found that the 
shape of the formation spectra is simple and seems to be 
suited to extract the binding energies and widths of the 
meson bound state because of the simple a — d cluster 
structure of ^Li. The size of the formation cross section 
can be more than half of the elementary cross section at 
the recoilless kinematics. 

This theoretical model is so simple that we can apply 
it easily to evaluate the mesic nucleus formation rate of 
other two nucleon pick-up reactions such as (tt, d) and so 
on. To do this, we should simply replace the elementary 
cross sections to those of the appropriate processes of the 
meson production like n + d ^ d+ heavy meson. In this 
model, however, the target nucleus is required to have 
large component of the quasi deuteron structure. Thus, 



we have considered the ^Li target as an example in this 
article. The calculated spectra shape are expected to 
have simple structure generally because of the existence 
of the quasi-deuteron in target nucleus and are suited to 
extract meson properties from the reaction spectra. 

In general cases, we should not assume the existence 
of the quasi deuteron in the target nuclei ^2d\ and we 
need to evaluate the emissions of the deuteron composed 
of two nucleons which are in the different single particle 
levels in the target. The deuteron can be formed from any 
pairs of proton and neutron in the target by the reaction 
of the meson production. In this process, however, the 
spectrum shape could be so complicated that it is difficult 
to extract meson properties. 

As for the actual experimental observations, the calcu- 
lated cross sections could be too small to find peak struc- 
tures in the inclusive missing mass spectra due to the size 
of the elementary cross section, and the coincident mea- 
surement detecting the particle pair emissions from me- 
son absorption in nucleus may be necessary to reduce the 
background. So far, the formation 77 mesic nucleus in the 
two nucleon pick-up (p,^IIe) reaction for ^^Al target was 
reported in Ref 15] by COSY-GEM. They performed the 
coincidence measurement with pn" pair emissions from 
rjN in nucleus. They reported 0.5 [nb] for the upper limit 
of the signal of the 77 mesic nuclear formation cross sec- 
tion [1^ at the energy where the elementary cross section 
is 77 [nb/sr] [2lj. In our example considered in this arti- 
cle, the effective number N^ff for the formation of meson 
bound state is about 0.5 for both shallow and deep po- 
tential cases around p^ — 1.6 [GeV/c] {Tp ~ 2.7 [GeV]) 
for the largest contributions as shown in Fig. [S] In this 
sense, we think the present result also have relevance as 
a guide for the actual experiment. 

We believe that it is quite important to find new re- 
actions suited to form the heavy meson bound states in 
nucleus to explore the various aspects of the strong inter- 
action symmetries at finite density by the mesic nuclei. 
In this context, the two nucleon pick-up reactions studied 
in this article are quite interesting, since we can satisfy 
the recoilless condition in this reaction for the formation 
of meson heavier than nucleon. 
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